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Abstract 

We compute analytically the diagonal quark number susceptibilities for a quark- 
gluon plasma at finite temperature and zero chemical potential, and compare with 
recent lattice results. The calculation uses the approximately self-consistent resum- 
mation of hard thermal and dense loops that we have developed previously. For 
temperatures between 1.5 to 5Tc , our results follow the same trend as the lattice 
data, but exceed them in magnitude by about 5 — 10%. We also compute the low- 
est order contribution, of order log(l/as), to the off-diagonal susceptibility. This 
contribution, which is not a part of our self-consistent calculation, is numerically 
small, but not small enough to be compatible with a recent lattice simulation. 



1 Introduction 

A lot of effort is presently devoted to understanding the properties of hot and 
dense matter from Quantum Chromodynamics. This is motivated in part by 
the ongoing experimental program on ultrarelativistic heavy ion collisions, and 
also by the progress in lattice gauge calculations which provide so far the best 
theoretical tool at our disposal to calculate from first principles the properties 
of the quark-gluon plasma. Recently however, it has been shown that results 
of such calculations could be remarkably well reproduced by weak coupling 
techniques when the temperature is larger than 2 to 3 times the transition 
temperature [1-3]. The purpose of this paper is to apply these techniques to 
the calculation of quark-number susceptibilities which have recently received 
considerable attention. 
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These quantities are interesting in several respects. First of all, they are to 
date about the only quantities tliat can be calculated on the lattice and pro- 
vide information about finite density [4-9] . (Recall that lattice calculations are 
still limited to zero chemical potential; susceptibilities involve derivatives of 
the thermodynamic functions with respect to /i, and their limit as — > can 
be computed on the lattice.) Susceptibilities have also been discussed lately 
in the context of heavy ion collisions, as they can be related to measurable 
fluctuations in conserved quantities [10-12]. However, the main question ad- 
dressed here is a theoretical one, namely, whether the recent lattice results 
in Refs. [8,9] can be explained within resummed perturbation theory, that is, 
without invoking genuine non-perturbative contributions. 

The lattice results [8,9] for the diagonal susceptibility x (cf- eqs. (2)-(3) below) 
at temperatures between 1.5 and 5Tc show a slow approach of the ideal-gas 
result from below, with deviations of about 15%. But the weak coupling ex- 
pansion of X completely fails to reproduce this behaviour. In massless QCD 
at = 0, this expansion is presently known to order g'^log{l/g) [13,14] : 




Xo ^ 2N 8 U 

-if(f)''-^°(«*) 

(with Xo = NT'^/3 the ideal gas value and Ng ~ N'^ — 1). Leaving aside the 
still undetermined g"^ contribution, one finds that the perturbative results lie 
above the ideal gas values for all temperatures of interest, and decrease with 
increasing T. 

We may relate this failure to that encountered in the perturbative calculation 
of the pressure [15]. In both cases, the difficulty of perturbation theory has its 
origin in collective phenomena which develop at the scale gT, and which in 
a strict perturbative expansion provide large contributions starting at order 
g^. To cope with this, various resummation schemes have been proposed [1- 
3,16-18]. Here, we shall use the one developed in Refs. [1-3], which focuses 
on the physical picture of the quark-gluon plasma as a gas of quasiparticles 
with properties determined by the "hard thermal loops" (HTL) [19,20]. This 
approach has proven to be successful in describing the lattice data for the 
thermodynamics of QCD down to temperatures as low as 2.5Tc. 

First lattice measurements of the off-diagonal susceptibility Xud have also been 
reported in Ref. [9]. This quantity vanishes for the ideal gas, so it probes 
directly the interactions in the system. We show that, when = 0, it is of 
order g^ log(l/g'), and we compute this lowest-order contribution in both QCD 
and QED. 
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2 Diagonal susceptibilities 



Quark number susceptibilities are generally defined as: 



where i, j are flavor indices, is the quark number density, and P is the 
pressure. With all quarks massless and /li — (as appropriate for comparison 
with the lattice results), all diagonal and all off-diagonal elements become 
equal, and we write 



Xij 



/i=0 



= X for i = j , Xij _^ = X for i 7^ j. (3) 



/i=0 



We shall evaluate the diagonal susceptibility x within the resummation scheme 
developed in Refs. [1-3]. This is based on the following expression for the 
fermion number density J\f in terms of the dressed fermion propagators A± 
(cf. eqs. (4.12) and (4.19) of Ref. [3]): 

AA=-47V /(04^{l^log^;' + Imlog(-A:i) + 

- ImE+ReA+ + ImE_ReA_}, (4) 

where A^^ = — [cj =l= (A; + E±)], E± are the corresponding sclf-cncrgics, and the 
plus (minus) subscript applies to fermions whose chirality is equal (opposite) 
to their helicity. The fermion self-energies and propagators are diagonal in 
flavor indices, and eq. (4) applies to each quark flavor i separately, but flavor 
indices are kept implicit. 

As in Refs. [1-3], we shall consider two successive approximations to the self- 
energies E±. The first is the HTL approximation where [19]: 



and is the plasma frequency for fermions, i.e., the frequency of long- 
wavelength [k — >■ 0) fermionic excitations {Cf = {N"^ — 1)/2N): 
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In this approximation, there is no mixing between quarks of different flavors, 
so the corresponding susceptibihties are diagonal even for /i 7^ 0. 

The resulting expression of the number density, denoted by Mhtl, is the sum 
of two contributions: Mhtl = -^htl '^■^htl^ where J^htl the contribution 
of the quasiparticle polesQ uj± = ±[k + T,±{uj±, k)], and J^htl is that of the 
Landau damping cuts at —k < uj < k. We have: 



^Ttl = N /^^{Tlog(l + e-[--(^)--l/^) 

^ 

1 + e-k-(fc)-/^]/r 

^T\og ^^^„(,_,v^ + (/i - -/i)}, (7) 

where the derivative is to be applied to the explicit yU dependence only, and 
not to that implicit in the dispersion laws and uj-{k) of the quasiparti- 

cles. Similarly, 



■'^ HTL 



00 , 2 k 



-N 



k'dk 



duj 



+ 



9/i 9/i 

u u 

X |arg[A; — uo + k)] — Im f^+{uj, k)'Re[k — 00 + k)]"^ 

+ axg[k + uj + t_{u,k)] - lmt_{uj,k)Re[k + uj + t^{u,k)Y^y 



The HTL approximation [19] contains the perturbative contributions of order 
g'^. This comes exclusively from the hard {k ~ T), "normal", branch uj+ and 
its asymptotic thermal mass = 2fcS_|_(ci; = k) = 2M^ : 

A/'^^^ = -T^/^M^. (9) 



However, there is no contribution in Mhtl- Such a contribution, denoted 
as A/'*-^\ comes entirely from the next-to- leading (NLO) correction 5M'^{k) = 
2fc Re (cj = k) to the asymptotic mass of the hard fermion [2,3]: 



where the NLO self-energy is given by the diagrams in Fig. 1. 

-"^ Charge conjugation (S+(a;, k) = w, A;)) exchanges the poles of A+ and A_ : 
A+ has two poles, one at positive uj, with energy and another one at negative 
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— > — > > > — 

Fig. 1. NLO contributions to ST, at hard momentum. Thick dashed and wiggly 
Hues with a blob represent HTL-resummed longitudinal and transverse gauge boson 
propagators, respectively. 

In contrast to the lowest order asymptotic mass M^, the correction 6M^{k) 
is a nontrivial function of the momentum [3] which can be evaluated only 
numerically. However, this function contributes to eq. (10) only in an averaged 
form [2,3] : 

^ jdkkduikW^ ^-2^' ^^^^^ ' ^''^ 



where rho, the Debye mass, is 

2rp2 

"6~ ^ 27r2 



ml = {2N + Nff-^ + ^^^„l (12) 



Thus, at a strictly perturbative level, it would be possible to reproduce the 
perturbative result for Af to order by replacing 2M^ = — > + 5M^ 
in eqs. (7,8) for the HTL approximation to jV. 

However, the correction 5M^ is negative, and for > 1 it is of the same order 
of magnitude or larger than the lowest-order asymptotic mass, apparently 
leading to a tachyonic thermal mass. As wc have argued previously [1~3], this 
problem is not specific to QCD, but can be studied already in simple scalar 
g'^ip'^ theory. There the perturbative thermal mass to NLO is = gf^T^(l — 
Sg/n). The corresponding one- loop gap equation, on the other hand, gives a 
monotonic function of g, which is well approximated by the quadratic equation 
[3] = g'^T'^ — SmT/n. Also a simple Pade resummation [1] = g^T'^/{l + 
Sg/n) gives reasonable approximations even for g ^ 1. In the following, we 
shall consider both prescriptions for including (11) into the asymptotic thermal 
mass, referring to them by "NLQ" and "NLP" , respectively. 



Lo, with energy —LO-{k); these go over to ±M as A: ^ 0. Correspondingly, A_ has 
poles at LO- and —0;+. See, e.g., Ref. [19]. 
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x/Xo 




0.1 0.2 0.3 0.4 0.5 

as 



Fig. 2. Quark number susceptibility normalized to its free-field value in SU(3) 
with Nf = 0,2,3 quark flavors as a function of Ug. Dashed lines give the per- 

turbative results through order (marked "(2)") and order as (marked "(3)"). 
The full line gives our result using HTL propagators, the dotted one a sim- 
pler quasiparticle model with momentum-independent mass M = M^o (both of 
which are A'^^ independent). Including next-to- leading order corrections to the 
asymptotic fermion mass through a quadratic gap equation gives the lines marked 
"NLQ"; using instead a simple Fade approximant gives the lines marked "NLF". 

3 Numerical evaluation 

The results of a numerical evaluation of x/Xo are given in Fig. 2 as a function 
of as- 

The HTL approximation gives results which are above those of first-order per- 
turbation theory (i.e., order g'^). Since the former does not include anything of 
the plasmon effect oc g^, the visible deviation is to be attributed to higher or- 
der contributions. A numerical analysis reveals that most of the enhancement 
is due to terms of order g'^, which also involve a logarithm : 

x2kU=o = N (^0.0431 ... X log ^ + 0.0028 • • •) (13) 

The coefficient of the log is by a factor of ~ — 0.52(iV^ — l)/iV^ different from 
that of the perturbative result (1). The correct coefficient will be restored by 
O [g'^ \og{l / g)T'^) corrections to (not considered here)[3]. 

^ The constant behind the logarithm (which is still unknown in perturbation the- 
ory) receives three-loop contributions which are beyond the ^-derivable two-loop 
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It is instructive at this stage to compare with the susccptibiUty of an ideal 
gas of massive fermions with mass equal to the asymptotic HTL mass, and 
which therefore contains the correct contribution of order g^. This reads (with 



Xr 



^ ' ° 



In contrast to eq. (13), this does not involve any logarithmic term at order g"^ 
,..U.;V^^.0.0216;V^. (15) 



Numerically, however, (14) happens to be rather close to the HTL expression, 
as can be seen from the dotted line in Fig. 2. 

At any rate, these order g'^ effects in either (13) or (14) are quite small com- 
pared to the more decisive order (7^-contribution. As we have seen, in the 
self-consistent density the effect of order g^ comes exclusively from the NLO 
correction to the asymptotic thermal mass. This introduces a (weak) depen- 
dence upon Nf, via the the Debye mass (12). As an estimate of this effect, we 
include it in the averaged form (11), for simphcity by a rescaling of M for all 
momenta. In order to get an idea of the theoretical uncertainties, we do so al- 
ternatively through a quadratic gap equation (NLQ) or through a (2,1)-Pade 
approximant (NLP). The corresponding numerical results for Nf = 0,2,3 
are shown in Fig. 2 by the various dash-dotted lines, with the formal limit 
Nf — corresponding to the quenched approximation of lattice gauge theory. 
As manifest on this figure, the inclusion of the order-^f^ contribution in our 
self-consistent calculation has a significant effect, although not as dramatic as 
in conventional perturbation theory. 

In Figs. 3 and 4 these numerical results are translated into plots of x/xo as a 
function of T/Tc using the recent determination of Tc/Aj-^ of Ref. [21] (which 
is found to differ significantly for quenched QCD and Nf = 2), together with 
a standard two- loop running coupling Q!s(/i). We vary the renormalization 
scale pL around p, — 27rT by a factor of 2. For an error estimate of the NL 
approximations, we in addition combine the (overlapping) results for NLP 
and NLQ. 

A completion of the g'^log{l/g) contributions, which is in principle possible 
within our approach and is left for future improvements, should decrease the 
NL results somewhat and presumably bring it nearer to the HTL result. 

approximation underlying the density expression (4). 
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A^^ = 0, fi = 'KT..A'KT 



T/T, 

Fig. 3. Comparison of our results for x/Xo in niassless QCD in the formal limit 
Nf = (using Tc/Ay^ =1.15 [21]) with the lattice results of Ref. [8] for quenched 
QCD obtained with quark mass m = 0.l2Tc on a lattice with Nf = 4 (no con- 
tinuum extrapolation). [The two rightmost lattice data are unpublished, private 
communication by S. Gupta.] 




Fig. 4. Comparison of our results for x/Xo in massless Nf = 2 QCD (using 
Tc/Ay^ = 0.49 [21]) with the lattice results of Ref. [9] obtained with quark mass 
m = O.lTc on a lattice with iVj = 4 (no continuum extrapolation). 

Also given in Figs. 3 and 4 are the recent lattice results of Refs. [8] and [9], 
respectively. These results involve finite but small quark masses, and, perhaps 
more importantly, are obtained for a lattice with only 4 sites in the temporal 
direction, and are still waiting for a proper continuum extrapolation. Our re- 
sults follow the same general trend as the lattice data (they slowly increase 
towards the ideal gas value), but exceed the latter by some +10%. (Remark- 
ably, this discrepancy is less pronounced for the physical case of dynamical 
fermions.) By contrast, the perturbative result to order g^, eq. (1), decreases as 
a function oiT/T^ in the range studied here and, actually, up to temperatures 
as high as T ~ lOOTc for Nf = 0, and even higher for Nf — 2. 
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4 OfF-diagonal susceptibilities 



The systematics of the diagrammatic contributions to susceptibilities (in par- 
ticular, the off-diagonal ones) can be clarified by referring to the symmetry un- 
der charge conjugation, or C-parity. Chemical potentials couple to the fermion 
fields in the same way as the A'^ component of an Abelian gauge field. Thus, 
when expanding a quark loop in powers of fi, one may attribute to each factor 
of fi the C-parity of the photon field, i.e., C = —1. Gluons attached to a quark 
loop in a colour symmetric state behave under permutations in the same way 
as photon insertions, and thus can be ascribed C = — 1 as well. C-parity con- 
servation forbids a photon to decay into two gluons: a colourless 2-gluon state 
is necessarily colour symmetric, and therefore C-even. However, a photon can 
decay into three gluons which are in a colour symmetric state, or in two gluons 
and an arbitrary odd number of photons, etc. In terms of chemical potentials, 
this means that a quark loop with two gluon external lines is necessarily even 
in n, while a quark loop with three gluon legs may generate also a term linear 
in n, which is then symmetric in the colour indices. 

The first perturbative contributions to the nondiagonal susceptibility x require 
two fermion loops connected by gluon lines. The diagram with just one gluon 
exchange vanishes by colour neutrality. The one with two gluon exchange 
is non-zero, but because the fermion loops are then even functions of n, it 
contributes to x only when yU 7^ 0, starting at order g^. One easily finds 

In fact, this is the same as Xij = aATf V% ^ j) with 

given m 

eq. (10). From the perspective of eq. (10), the mixing between different flavors 
is induced by the resummation of quark loops along the soft, internal, gluon 
lines in the diagrams in Fig. 1. 

However, when all chemical potentials vanish, the lowest-order diagram con- 
tributing to X is the " Ibugblatterf ' [22] diagram shown in Fig. 5a. 



This diagram is superficially of order g^, but when calculated with bare gluon 
propagators, it develops a logarithmic infrared divergence in the electrostatic 
sector, because, for static external gluons, the quark loop induces an effective 
local vertex corresponding to 

^^TrAlY^£^P^{mf,fi,,T) (17) 
where Pq is the ideal gas pressure of a fermion with mass nii and chemical 
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(a) (b) 

Fig. 5. (a) Lowest-order diagram in the thermodynamic potential that contributes 
to off-diagonal susceptibilities dMijdiij at /i = 0. (b) Corresponding diagram in the 
effective theory for the electrostatic modes. 

potential For massless fermions this reduces to [23-25] : 

I^IV AlY^^x,^ r^^d'^'-'AlAlAl Y: l^i (18) 



{A^ = Ap"', TrtH'' = 5"^/ 2). We expect Debye screening to cut off this 
divergence at the scale rhjj ^ gT, with the upper scale in the logarithm, of 
order T, set by the thermal distribution. This gives a contribution to x of 
order log(l/5'), which is the leading order effect for g small enough. We now 
compute its coefficient. 

In the imaginary time formalism, the infrared divergence is isolated in the 
static Matsubara sector. The original diagram in Fig. 5a is then identified with 
the two-loop diagram in Fig. 5b. Formally, this is the second order perturbative 
correction /2 to the free energy f — — log Z of a 3-dimensional scalar field with 
effective action 

Se = J d'x^A^,{-V' + ml)A^, + ^Y,|^J^-^d'^''A^,A',A^„ (19) 



where ^^(x) = y/TA^{uj^ = 0,x), ttid is the Debye mass (12), and the in- 
teraction term is now purely imaginary, as a consequence of the continuation 
A^ — > IAq to imaginary time. Denoting by Sj the interaction term in eq. (19), 
we have /2 = —{S]/2)o, which is positive, a consequence of the interaction term 
being purely imaginary. A direct calculation yields: 



= M'^'^d^"^- (^E/^.^j /-^I^oo(A;)i^oo(g)/^oo(|k + q|),(20) 

where Dqo = l/(/^^ + ml) and d'^'^'d'''"' = {N^ - 1){N^ - A)/N. The above 
integral has a spurious ultraviolet divergence, which comes from the restriction 
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to the static Matsubara modes, and, in the absence of the Debye mass, it would 
be also divergent in the infrared. This yields 



/ 



^d?k (Pq 



1 1 A 1 , 1 

Idtt^ niD Idtt"^ g 

where the upper cut-off A eventually gets replaced by T upon inclusion of the 
nonstatic Matsubara modes. 

Putting everything together and returning to the original 4-dimensional gauge 
theory at finite temperature, the above estimate for /2 translates into the 
following, negative, contribution to the pressure: 



Prom eq. (22) we finally deduce the leading-order term for nondiagonal x '■ 
^ ^ ^ ' ^ Mog-. (23) 



This logarithmically enhanced contribution vanishes in SU(2) gauge theories — 
though not in QED. The (ultrarelativistic) QED result is obtained by replacing 
^abc^abc _^ iQ eq. (20) (cf. eq. (18)), and xo ^ ^73, yielding 



X/Xc 



QED 



87r6 



log(l/e) ~ -4(Qi/7r)Mog(l/Q;) 



(24) 



In (massless) QCD, the leading-order contribution to x is 

^ ~ -— f^Vlogi ~ -^aMog(l/o;,). (25) 



While in QED it is plausible that the leading-log term with its negative coeffi- 
cient dominates over other contributions oc a^, this is more uncertain in QCD, 
where log(l/Q;s) is much smaller. But assuming that the unknown constant 
behind the log is roughly of order 1, the off-diagonal susceptibilities appear to 
be rather tiny in QCD (although they would tend to become more important 
for larger N). In fact, a most recent lattice study of nondiagonal susceptibili- 
ties in Nf — 2 QCD [9] has found only values consistent with zero, but within 
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statistical errors that are < 10 ^, whereas the natural order of magnitude of 
(25) is given by ^[a,(27rr)]3 ^ 10"^ for T - 

The lattice calculations in Ref. [9] have been performed with finite quark 
masses down to m/Tc — 0.1. However, this should not lead to any noticeable 
reduction, because the first m/T correction in (17) is only quartic, leading 
to a correction factor ?a (1 — 0.06187m"^/T^) in the final result (25), so the 
extreme smallness of the lattice result of Ref. [9] remains a mystery for now. 



5 Conclusions 



To summarize, we have presented an analytical calculation of the diagonal 
quark number susceptibility in hot QCD within an approximately self-consis- 
tent resummation of perturbation theory. Our (non-perturbative) formulae 
include completely the pcrturbative contributions of 0{g^) and 0{g'''). For 
temperatures between 1.5 to 5Tc, our results show the same general trend as 
seen on the lattice - namely, a slow increase towards the ideal gas results from 
below -, but with absolute values which are slightly, but systematically, above 
the lattice data, by 5-10% in the case oi Nf — 2. This deviation is somewhat 
larger than that of our analogous calculations of the entropy density [1-3]. 
However, given that a continuum extrapolation of the lattice data for quark 
susceptibilities is still missing, it remains to be seen whether there are sizable 
higher-order pcrturbative contributions not captured by our approach or even 
important non-perturbative phenomena, as speculated in Ref. [9]. 

We have further computed the off-diagonal susceptibility to lowest non-trivial 
order in perturbation theory, that is, to order g^log{l/g). The result turns 
out to be remarkably small, but not so small, however, to explain the corre- 
sponding lattice result of Ref. [9] , which, surprisingly, is consistent with zero 
with statistical errors < 10~^. This discrepancy certainly calls for more inves- 
tigations and more lattice data. 
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